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Abstract 

The time evolution of systems relaxing towards thermal equilibrium is exam- 
ined near the critical temperature T c , with special attention paid to the role 
of the initial value m; of the order parameter <fi. To this end, the n-component 
model A for a cube of length L is investigated. The common belief that all 
memory of m; is necessarily lost after a microscopic time span is shown to be 
unfounded. General arguments and the exact solution of the limit n — > oo 
show that mi leaves its traces in both the linear and nonlinear long-time re- 
laxation of (f) near or at T c . Specifically for linear relaxation near T c , or at T c 
with L < oo, the amplitude of the exponential decay depends on m; and the 
short-time exponent 9' = (x{ — x^/z, provided fj ~ m i z ^ x ' is comparable to 
or larger than other time scales. Here X{ is the scaling dimension of m;, z is 
the dynamic bulk exponent, and is the usual equilibrium scaling dimension 
of </>. 

PACS numbers: 05.70.Jk, 05.70.Ln, 64.60.Ht, 75.40.Gb 
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How does a thermodynamic property of a system at temperature T above, but near 
its critical point T c , relax to zero with time t? This is a long-studied problem of critical 
dynamics, which we wish to re-examine here. The main question we wish to answer is 
whether and to what extent the initial condition affects such relaxation processes on long 
time scales. 

The stochastic models used in the theory of dynamic critical phenomena [1] to represent 
the various dynamic universality classes, by construction, possess a number of properties 
such as detailed balance which guarantee that the system relaxes from almost any initial 
configuration towards the thermal equilibrium state. The common belief is that all memory 
of the initial state gets lost on a microscopic time scale. As we shall show, this is not in 
general the case: Under appropriate (realizable) conditions (to be discussed below), the 
initial condition may well leave its fingerprints on the relaxational behavior on long time 
scales. 

To become more specific, let us consider the following relaxation process: A ferromag- 
netic system with non-conserved n-component order parameter <p and purely relaxational 
dynamics (model A of Ref. [1]) is prepared in a macrostate with homogeneous magnetiza- 
tion rrii at temperature To 3> T c and magnetic field h > 0. At t = the magnet is rapidly 
quenched to a temperature T > T c and h is switched off. Then the system is let to evolve. 

Immediately after the quench, up to some microscopic time scale t m - lc , the time evolution 
clearly should be non-universal, depending on both the microscopic details of the interactions 
in the sample and the initial condition. For t — Y oo, the same unique equilibrium state is 
approached, irrespective of the precise initial condition. This suggests that for sufficiently 
large t 3> t m - lc the process should enter a regime of universal behavior, involving only a few 
macroscopic time scales. 

One well-known such scale is the characteristic time t T ~ t~ vz , where r = (T — T c )/T c ; 
in our case it may be defined as the linear relaxation time [2] of the order parameter 
m(t) = (</>(x,f)) in a bulk system. For a long time it had been believed that only two 
time regimes with universal time dependence could be distinguished [1]: a regime t 3> t T of 
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linear relaxation, m{t) ~ e~ i//tr , and a regime t m ; c <C i <C i T of nonlinear relaxation [3,4] with 
m(t) ~ t~^' vz . The problem with this simple picture is that it precludes the possibility of 
any initial-time effects on macroscopic time scales. This is unjustified since the initial value 
rrii yields an independent time scale that may become arbitrarily large for sufficiently 
small rrii. In fact, if we assume that m; scales with a dimension Xi under dilatations, then 
ti ~ m i z l x \ This suggests that the above simple picture is incomplete and needs to be 
corrected. 

A crucial step in this direction was made by Janssen et al. [5,6] who showed that, contrary 
to previous assumptions, X{ is a genuine new exponent, generally different from the scaling 
dimension x^ of <f>(t > 0). As a consequence, m{t) was found to exhibit a non-trivial short- 
time behavior of the form m(t) ~ t 6 ' for t m - lc < i < < i T] with 8' = (xi — x^)jz. For 
model A below its upper critical dimension d* = 4, Xi turned out to be larger than x^, so 
that m(t) actually increases in this regime. 

Ref. [5] and subsequent work [7] mainly dealt with the critical case (t T = oo). To 
complement these findings we here focus our attention on regimes in which t ^> t m - lc is large 
compared to the linear relaxation time. While our main interest is in the bulk case, we will 
also discuss the finite-size case of a cube of length L with periodic boundary conditions, so 
that the special case t T = oo, L < oo, treated in Ref. [8] is included. For this special case we 
have shown there that the initial condition also affects the linear relaxational behavior on 
long time scales. Using general scaling arguments, we show below that this phenomenon is 
not restricted to finite-size systems at T c but can be observed quite generally in the case of 
linear relaxation. Our main findings will be corroborated by exact solutions for the n — > oo 
limit of the mentioned n-component model A with L < oo. 

To develop our argument we begin by recalling the scaling form [8] 

m(t, r, L, rrn) « C m t~^ vz T(t/t T ,t/t L ,t/t t ) (1) 

phenomenological scaling considerations lead one to anticipate on large length and time 
scales. Here C m is a non-universal amplitude. Further non-universal metric factors appear 
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in t T , ti, and t{. The function T is universal. 

In Ref. [8], the form of J 7 (0 , $ l , i n various limits was discussed. For example, the 
anomalous short-time behavior ~ t 6 ' could be traced back to the fact that .F(0, 0, $) exists 
and behaves as 

^(0,0, tf) ^T ha ^ z (2) 

as d = tjti — y [9], where Tbd is a universal constant. 

In the long-time limit, with which we are concerned here, m(t) must decay exponentially, 
unless t T = II = oo. That is, as t — > oo with fixed t T , t^, and f;, it should vary as 

m{t,t T ,t L ,t i )lC m ziM 00 e- tltE . (3) 

Compatibility with (2) requires that the amplitude M m and the linear relaxation time 
have appropriate scaling forms. The former in general depends on r, L, and rrii] it can be 
written as 

- /../».: = t-^ V2 A(t L /t T ,ti/t T ) . (4) 

On the other hand, t% should be independent [11] of and satisfy the relation 

t E (r,L) =t T T(t L /t T ) . (5) 

Since t T = ^(t, oo) and t E = t^(0, L), the function T(x) must have the following asymptotic 
properties 

T(x — )■ oo) ?s 1 and T(rr — >• 0) « rr . (6) 

Turning to the universal amplitude function „4(:c,y), we first note that the bulk limit 
x — y oo of (4), and hence 

.4(:e -)• oo,*/) « A(y) , (7) 

should exist. At criticality, (4) should become 
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« t^'^EXUlH) ■ (8) 

This implies the behavior 

A(x,w/x) « x-W vz B c {w) (9) 

as a; — y with fixed w = ti/tL. If both t T and ^ are large compared to ti, it seems reasonable 
to expect (in complete accordance with our exact results below) that any dependence on 
rrii should drop out of the limiting form of m(t) and hence of M^. From this we infer the 
asymptotic properties 

A(x,0) = AM and B c (0) = B coc , (10) 

in which B coc ^ is a finite universal number. 

Of particular interest are the cases with t{ large compared to t T , to ti, or to both. For 
those our exact results below and in Ref. [8] reveal a dependence of on wi;, suggesting 
that this dependence should be linear. In order that this behavior complies with (4) we 
must have 

A(x,y -}> oo) « y-^ z Ai(x) . (11) 

From this we immediately deduce the limiting form 

Moo ~ const m ; tf Ai(t L /t T ) (12) 

for t{ — y oo with fixed t T and t^. This is a central result of this Letter. It shows that the 
amplitude depends indeed on the initial value rrii, and that the short-time exponent 6' 
governs the r -dependence of M m . For the function A{ we anticipate the limiting properties 

Ai(x — y oo) « „4 ife and ^(rr — 5- 0) « „4 ic . (13) 

The first ensures that both above features survive the bulk limit L — y oo; the second implies 
that (12) reduces for t T 3> ti 3> t^ correctly to the r = form 
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Moo ~ const nii t e ^ Ai c (14) 

discovered in Ref. [8]. 

The scaling forms predicted above may all be verified explicitly in the limit n — ¥ oo of 
our model. This is defined by the Langevin equation 

A" 1 a t 0(x, t) = (- A + r + r c + |- <fA<f> + C, (15) 
V on J 

in which r c is the usual static critical value of the bare squared mass, while ( is a Gaussian 
random force with zero mean and variance 

<C«(x, *) (p (x', t')} = 2 8 Qf) 5(x - x') 5(* - *') (16) 

We assume that 2 < d < 4, so that an ordered equilibrium state exists and hyperscaling is 
valid. In the limit n — > oo the model becomes Gaussian with a time-dependent susceptibility 
x(t)] its dynamics is completely described by the self-consistent set of equations 

X{t)d t m{t) = -Xm{t) (17) 

and 

x{ t)-i = T + T c + ^{C(t) + m(tf] (18) 

with 

C(t) =]>]C(q;M) , (19) 
q 

where m(i) now means the rescaled magnetization n _1//2 |(0(x, t))\ with initial value m(t = 
0) = wij. For brevity all arguments of m(t,r, L,nii) but t are omitted. The sum in (19) 
extends over discrete momenta q = 2 7rm/L with m £ The Fourier-transformed auto- 
correlation function C(q;i,i) may be expressed in terms of the response propagator 

G(q;M') =8{t-H) exp (-A jf dt" [q 2 + x(0 _1 ]) (20) 

as 
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roc 

C(q;t,t)=2X dt' G(q; t, t') G(q; t, t') . (21) 
Jo 

When (17) and (18) are combined, %(t) can be eliminated to obtain a linear integro- 
difFerential equation for f(t) = l/m 2 (t) that can be solved by Laplace transformation [12]. 
The result for the Laplace transform f(s) = / °° dte~ st f(t) can be written down in explicit 
form. Ignoring corrections to scaling, we find 

f(s) ~ (2XA^L d EillM+lhlllL (22) 
/ (s) - (2A A d ) L h (sL2 /2A) _ ^ L d-i T j Ad ( 22 J 

where 



with 



oo 

,2 , 



HO = £ <r n /4? - 1 



(24) 



A d = -(ATi)- d l 2 Y(l-dl2) , and D g = 6/g. 

The result can be cast in the expected scaling form 

f{s) = const t^^Gist^tJtL^/tL) , (25) 

with 1 + 2(3 /vz = d/2 for the n = oo case considered here. The relaxation times t T and 
[normalized in accordance with (3) and (5)] are explicitly given by 

1 (D„t\- vz L z 

with the n = oo values vz = 2/(d — 2) and z = 2. The parameter xq denotes the zero of 
h(x). Upon taking ti as 

D„ yd 2) —zlxi lr\n\ 

t t = — l m . (27) 
with the n = oo value z/xi = 2, one can derive from (25) the scaling function 



r( \ 2w/(d-2) + l/u 

y(u.v.w) = — — . ... , . 28 

v ' ' ' h(x Q u/2) - (x /vY d - 2 )/ 2 v ' 

For general values of t T , II, and ti the Laplace back transform of (22) — and hence the 
scaling function T — cannot be calculated analytically but must be computed numerically 
[8]. However, TbJyff) = ^(O, 0, the scaling function pertaining the bulk critical case, can 
be obtained in closed form. It reads 

™ = {jTi) ^ (29) 

For d < 1, T hc « tf 1/2 , in accordance with (2). In the large-?? limit, T hc = 1 - l/2t? + 0(t?" 2 ) , 
i.e., the leading term is independent of d. This bears out that the resulting asymptotic long- 
time form of m(t), the well-known nonlinear relaxation, becomes indeed independent of the 
initial value m;. 

The asymptotic long-time behavior (3) of m{t) in the case t T + > 0, and the scal- 
ing function A of (4), can also be obtained analytically, employing standard results about 
Laplace transformation [13]. The asymptotic exponential decay follows from the contribu- 
tion of the pole of f(s) with largest positive real part, which in turn may be identified with 
2/tg. One thus finds that the scaling function T(x = tL/t T ) of (5) is the solution to 

h(x Q x/T) - (x Q xf- 2 ^ 2 = . (30) 

Noting that h(xo) = and h{x -¥ oo) ~ rr^ 2-1 , one easily verifies that the limiting properties 
(6) are satisfied. 

The result for the scaling function A(x,y) reads 

= { ™ { $vS <3i) 

with the n = oo exponent value 6' = (4 — d) /4 and where we have chosen the normalization 
such that A(x — > oo,y — > 0) — > 1. Our result bears out the fact that the amplitude 
defined in (3), in general, depends on t{ and hence on m;. Memory of the initial value is lost 
only by decay. 



Starting from (31) all the asymptotic cases discussed above can be studied in detail. 
Especially the scaling functions Ab, B c , Aoc, and Ai introduced in (7), (9), (10), and (11), 
respectively, can be calculated explicitly. For y = in (31), we obtain 



Aocix) = (X X) 



2h'(x x/T) 



1/2 



(32) 



(d - 2) r 

In this case the dependence on the initial field drops out, only subleading terms in the 
asymptotic expansion depend on the initial conditions. 

In the complementary limit, y — > oo, we identify by comparing with (11) 

Ai(x) « 2- 1 / 2 (x Q x) 6 ' [h'(x x/r)] 1/2 , (33) 

and the limiting forms (12) and (14) can be verified. Hence, in the regime t{ 3> t T , the 
amplitude M m indeed depends linearly on m;. Further, the r or Z-dependence of for 
ti th^> tr or ti t T ^> II, respectively, is governed by the universal short-time exponent 
8' . For the latter case it is essential that 8' also occurs in the scaling function A(x, y). 

So far we have focused our attention on the magnetization, the one-point function of the 
field theory. A natural question to ask is whether the above long-time dependence on m; 
also occurs in other quantities, in particular, in two-point functions. For our model with 
n = oo the result for the response propagator reads 

G( q] t,t')=8(t-t')^e-^^ . (34) 

Thus, in the t,t' — ¥ oo limit it is given by 

G(q; t, t') « 8{t - t') e -(M 2 +i/t E )(t-t') ^ ^ 

that is all dependence on m; drops out. A similar result holds for the two-point cumulant. 

A second obvious question is whether similar long-time traces will also occur for other 
models. Since our scaling arguments have been fairly general, we expect this to happen 
for those universality classes that have a non-trivial critical short-time behavior — i.e., for 
models with a non- vanishing short-time exponent 8'. As shown in Ref. [7], interesting critical 
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short-time behavior also occurs when (i) the order parameter is coupled to a conserved 
density (model C) and (ii) when one has reversible mode coupling (models G and E). It 
would certainly be interesting to extend our above analysis to these cases. 

In summary, we have shown that, against the common belief, initial conditions affect 
critical relaxation processes even in the long-time regime. As expressed in (4) and explicitly 
verified for the large-n limit, the amplitude of the linear decay is governed by the initial 
value of the order-parameter field as long as f;, the initial time scale, is comparable or larger 
than other time scales. As critical dynamics is now in reach of computer simulations [14], it 
should be also possible to verify our results numerically. 
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